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We present an improved uvby-metallicity relation calibrated for F,G, and 
early K dwarfs, and an analogous uvby-T e /j relation, both derived using a 
Levenberg-Marquardt minimization scheme. Our calibrations are based on 1533 
stars which appear in both the Cayrel de Strobel (2001) metallicity compila- 
tion, and in the Hauck-Mermilliod (1998) catalog of uvby photometry. We also 
examine the speculative possibility of using uvby photometry to produce a uvby- 



O ■ planeticity calibration. We conclude that while there is likely no strong photo- 

^ ■ metric indicator of the presence or absence of short-period planets, the possibility 

^ ■ of a spectroscopic indicator of planeticity is well worth examining. 

CLi 1 Subject headings: stars: planetary systems — stars: metallicity 
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1. Introduction 

> : 

^ . The number of known extrasolar planets is approaching 100, and more planets are being 

discovered every month. 1 A remarkable empirical correlation has emerged from the aggregate 
of planets: The parent stars of the detected extrasolar planets appear to be significantly 
metal-rich in comparison to field FGK dwarfs in the solar neighborhood. This correlation 
was first extensively discussed by Gonzalez (1997), and has been further studied by a number 
of authors; a partial list might include Gimenez (2000), Laughlin (2000), Santos, Israelian & 
Mayor (2001), Murray & Chaboyer (2001), and Reid (2002). These studies all agree that the 
fraction of stars containing readily detectable extrasolar planets (e.g. having K > 15 m/s, 
and P < 5 yr) increases substantially with increasing stellar metallicity. A search strategy 
geared toward the most metal-rich stars should thus detect short-period planets at a far 
greater rate than a conventional volume-limited survey. 



For the latest catalog of planets, see http://www.exoplanets.org 



- 2 - 



One might argue that a metallicity-based planet search strategy will result in a haphaz- 
ard accumulation of short-period planets, while simultaneously skewing the overall census of 
planets in the solar neighborhood. Indeed, it is important to continue with volume-limited 
surveys in order to round out the unbiased statistical distribution of planetary systems. The 
discovery of additional short-period planets, however, has important ramifications. For the 
following reasons, it is desirable to locate as many short-period planets as possible: (1) 
The presence of a short-period planet is an excellent indicator that additional planets are 
present in the system. Currently, for example, 5 of 12 short-period planets observed at Lick 
Observatory show evidence for additional companions (Fischer et al 2001). Recently, it has 
been realized that some multiple planet systems such as GJ 876 allow for unambiguous de- 
termination of planet masses and orbital parameters (Laughlin & Chambers 2001, Rivera 
& Lissauer 2001). (2) The current census of short-period extrasolar planets shows an inter- 
esting concentration of objects in the period range between roughly 2.98 and 3.52 days (8 
out of 20 planets with periods P < 20 d). While this pile-up is probably a consequence of 
disk-protoplanet tidal migration (Lin, Bodenheimer & Richardson 1996), the stopping mech- 
anism is not at all well understood. It is thus worthwhile to find more short-period planets to 
better delineate the minimum planetary period. (3) Geometric arguments indicate that 10% 
of Hot Jupiter-type planets transit the face of the parent star. A large increase in the rate 
of short-period planet detections thus translates into a significant increase in the number of 
transits discovered. Transiting planets permit accurate mass and radius determinations, and 
are extremely important data points within the overall theory of giant planets. 

High-resolution spectroscopy provides the most precise method for flagging the metal- 
rich stars which have a high probability of bearing a detectable planet. Unfortunately, 
the number of field stars with spectroscopically determined metallicities is relatively small. 
The most recent compilation from the literature by Cayrel de Strobel, Soubiran, & Ralite 
(2001; hereafter CSR) contains 6534 metallicity determinations for 3356 stars. The majority 
of single, metal-rich, main sequence and subgiant FGK stars in the CSR compilation are 
already under radial velocity surveillance. A larger pool of high metallicity candidates is 
needed. 

A study by Laughlin (2000) used the Hipparcos database (Perryman et al 1997), the 
Hauck-Mermilliod (1998) compilation of uvby photometry, and the Schuster & Nissen (1989; 
hereafter SN) uvby-[Fe/H] calibration to to identify 200 metal rich stars from a parent 
population of roughly 10,000 K,G, and late F-type main sequence stars (having 7.8 < V < 10) 
in the immediate (d < lOOpc) solar neighborhood. Nineteen of these stars were added to 
the Keck Radial Velocity survey in the Summer of 2000. Three of these nineteen stars soon 
proved to be accompanied by planets (Vogt et al 2001, Butler et al 2002). Metallicity-based 
search strategies are therefore clearly optimizing the planet-detection efficiency 
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The SN calibrations for F and G spectral types were based on 219 calibration stars 
spanning a range of metallicities from —3.5 < [Fe/H] < 0.4. This very useful calibration, 
which has an overall standard deviation of 0.16 dex per star, is wholly empirical, and has 
the advantage of using only the standard uvby indices, rather than the differential indices, 
5mi = rn^ yades — raf ar , and 5c\ = c ^ yades _ d tar , which are used in the earlier calibrations of 
Crawford (1975) and Olsen (1984). The SN calibration has been widely used to investigate 
metallicity trends among solar neighborhood stars (e.g. Rocha-Pinto & Maciel (1996), or 
Haywood (2001). As recently emphasized by Twarog, Anthony-Twarog and Tanner (2002), 
however, the SN calibration systematically underestimates [Fe/H] by ~ 0.2 dex for stars with 
spectroscopic metallicities [Fe/H] > 0.1. The SN calibration is thus not optimally suited for 
compiling an accurately ranked catalog of high metallicity stars. 

Our main goal in this letter is to make the uvby-metallicity calibration more discerning 
and accurate among the highest metallicity stars. We also use this forum to briefly introduce 
the concept of a "planeticity calibration" , in which stellar properties are correlated directly 
with the detectable presence of an extrasolar planet. 



2. New Photometric Calibrations 



To produce our new calibration, we we queried the CSR database for all stars which 
have (1) Hipparcos parallaxes, (2) uvby measurements in the Hauck-Mermilliod (1998) com- 
pilations, (3) M v < 1.0, and (4) lie within 100 pc of the Sun (to reduce the problems posed 
by reddening). This subset of the CSR catalog contains 1533 measurements of 664 stars. 
Using a Levenberg-Marquardt scheme (see e.g. Press et al 1992), we then optimized the 
coefficients, A n , of the most general third-order polynomial in the Stromgren indices (b — y), 
mi, and c\ so as to best reproduce the spectroscopically determined [Fe/H] values. Because 
the Levenberg-Marquardt routine provides local convergence, it is prudent to start the itera- 
tive sequence with different sets of initial guesses for A ni to better insure that A n , the set of 
final coefficients produced by the routine, corresponds to the global minimum. We ran the 
Levenberg-Marquardt procedure using (1) A ni from the SN "F star" calibration (ignoring the 
logarithmic terms), (2) A ni from the SN "G star" calibration, and (3) by setting all A ni to 
unity. In each case, the routine converged to an identical set of final coefficients, suggesting 
that a global minimum has been achieved. Our calibration is: 



[Fe/R} phot = - 10.424602 + 31.059003(6 - y) + 42.184476mi + 15.351995ci 

- 11.239435(6 - yf - 29.218135m? - 11.457610c? - 138.92376(6 - y)m x 
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- 52.033290(6 - y)a + 11.259341miCi - 46.087731(6 - y) 3 + 26.065099m? 

- 1.1017830c? + 138.48588(6 - yfm 1 + 39.012001(6 - yfc x 
+ 23.225562m?(6 - y) - 69.146876m? Cl + 20.456093c?(6 - y) 

- 3.3302478c?mi + 70. 168761 (6 - y)m 1 c 1 

(1) 

Figure 1 shows our estimates of [Fe/H] pftot arising from this calibration plotted versus 
[Fe/H] spec for all 1533 calibration measurements. 

The CSR database has been compiled from a diverse array of metallicity determinations, 
some of which have reported uncertainties, and some of which do not. Furthermore, many 
of the 664 stars which survived our cuts for inclusion into the calibration are represented 
by multiple spectroscopic measurements. We circumvented the non-uniform property of 
the CSR compilation by assuming that every spectroscopic metallicity measurement has an 
associated uncertainty of 0.05 dex. With this assumption, each of the 1533 measurements 
contributes an equal weight to the calibration. We compute a reduced x 2 statistic of 6.68 for 
the overall fit, which is superior to the reduced x 2 value of 9.33 which results from comparing 
[Fe/H] p?i0t obtained with the SN calibrations with the spectroscopic metallicities. 

In Figure 2, we plot the distribution of metallicity differences between the photometric 
and spectroscopic determinations, [Fe/H] p /j Gt — [Fe/H] spec , for both our calibration and also 
for the SN calibration. Both distributions are well represented by Gaussians. The new 
calibration is superior to SN in that it has a narrower half-width at half maximum (0.10 
dex as opposed to 0.13 dex), and is centered more closely to zero (-0.027 dex as opposed 
to -0.049 dex). Figure 3 shows a similar diagram in which we plot only the results of 
measurements which yielded a spectroscopic metallicity [Fe/H] > 0.0. Here again, our 
calibration is significantly superior to SN, which is not surprising, as very few high metallicity 
stars were included in the SN calibration, which was geared primarily toward studies of 
population II stars. Our calibration still tends to underestimate the metallicity of the most 
metal-rich stars in comparison to spectroscopic determinations, but the discrepancy has been 
significantly reduced. As a specific example, the planet-bearing star HD 145675 (14 Her), 
which is the highest metallicity star in our subset of the CSR catalog, with [Fe/H] = 0.5 
(Gonzalez, Wallerstein & Saar 1999), receives a metallicity estimate of [Fe/H]=0.47 from our 
calibration, and [Fe/H] =0.13 from SN. 

Our metallicity calibration should provide a useful sieve for selecting candidate stars 
to be added to radial velocity surveys. We note here a handful of class IV and V stars for 
which (1) our new calibration finds [Fe/H] p h ot > 0.3, which (2) do not appear in the Laughlin 
(2000) catalog, and which (3) have no previously announced planet: HD 88176, HD 112164, 
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HD 15942, HD 144585. 

In addition to providing spectroscopic metallicities, the CSR catalog also lists effective 
temperature estimates for each of the cited entries. Following the same procedure and 
applying the same cuts as for our our [Fe/H] p ft ot calibration, we have also derived a uvby—T e ff 
calibration, 



[TefrW = 4434.8364 + 29839.391(6 - y) - 13967.518m! - 3041. 0698c! 

- 98624.604(6 - yf - 8963.3232m? + 10620.269c? + 63289.783(6 - y) mi 

- 10099.516(6 - y) Cl + 16204.483miCi + 91256.827(6 - y f - 45093.846m? 

- 10580.081c? - 122571.68(6 - yfm x + 16993.421(6 - yf Cl 
+ 111551.98m?(6 - y) - 43418.471m?ci - 5856.0002c?(6 - y) 
+ 15492.020c?mi - 127.04874(6 - y)m lCl 

(2) 



If we assume that the average measurement uncertainty of the spectroscopic T e ff deter- 
minations is 50K, the reduced x 2 of this calibration is 3.17, indicating that uvby photometry 
is quite adequate to providing effective temperature estimates for a given star. As before, 
we also construct a distribution of differences T e ff hot — T e ff spec , and construct a gaussian fit 
to the distribution. The half-width at half maximum for this distribution is 80.4 K, whereas 
the center lies at -11.5 K. 



3. Toward a Planeticity Calibration 

Given the demonstrated connection between the metallicity of a parent star and the de- 
tectable presence of a planetary companion, and given the power of the Levenberg-Marquardt 
minimization technique in producing empirical calibrations, one can ask a rather provacative 
question: Is it possible to construct a "planeticity" calibration in which the chance of finding 
a readily detectable planet around a given star is estimated directly from uvby photometry? 

Cumming, Butler and Marcy (1999), have published a detailed analysis of 11 years of 
precision radial velocity measurements for 73 nearby stars. Within this list of stars, there 
are currently 10 announced planets. We assign a planeticity V — 1.0 to stars with planets, 
and V = 0.0 to stars without. We then use the Levenberg-Marquardt method to produce 
a uvby-planeticity calibration in exact analogy to producing the metallicity and effective 
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temperature calibrations. This exercise leads to the following empirical calibration: 

[Vphot] = 168.92038 - 1367.1865(6 - y) + 1215.0485mi - 625.29367c! 

+ 2450.1524(6 - yf + 251.91349m? + 686.49412c? - 3024.1263(6 - y)mi 

+ 3776.9231(6 - y)a - 3394.6622m lC i - 838.40082(6 - yf - 347.90246m? 

- 238.02961c? + 698.73217(6 - yfm x - 3804.6842(6 - yfc x 

+ 802.37174m?(6 - y) - 454.17907m?ci - 2429.7350c?(6 - y) 

+ 2468.8122c?m! + 4574.5147(6 - y)m 1 c 1 

(3) 



The results of this calibration are shown in Figure 4, and indicate that the planeticity 
calibration is not particularly convincing. Nevertheless, it is instructive to test the relation 
on a larger set of stars. We queried the Hipparcos catalog for stars lying within 50 pc of 
the sun, retaining those for which (1) uvby photometry exists, (2) fall in the range 0.8M Q < 
M* < 1.5M Q (using the stellar evolution models of Allende Prieto & Lambert 1999), and 
which (3) were not among the 73 Lick Survey stars used to build the planeticity calibration 
given above. These cuts left 2056 stars, of which 46 have reported planets. We expect 
that the majority of the remaining 2010 are currently under radial velocity surveillance, 
and that relatively few easily detected planets will emerge from this aggregate. Using our 
calibration, we compute an average planeticity of 0.271 for the 46 known planet-bearing 
stars, in comparison to an average planeticity of 0.113 for the 2010 stars without reported 
planets. 

We conclude that uvby photometry likely holds little planeticity information beyond 
that provided by the well-recognized metallicity proxy. We are very optimistic, however, 
regarding the possibility of producing a planeticity calibration which is based on spectro- 
scopically determined differential metal abundances within a particular star. Recent work, 
in particular by Gonzalez et al (2001), indicates that planet-bearing stars may have system- 
atically unusual abundance ratios for [Al/Fe], [Mg/Fe], and [Na/Fe]. Debra Fischer (2002, 
personal communication) has recently completed detailed spectroscopic abundance analyses 
of more than 500 stars which are currently maintained in the Keck and Lick planet search 
lists. In a future project, we are hoping to replace the uvby indices with spectroscopically 
determined differential abundances, and thus derive a more definitive planeticity calibration. 
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Figure Captions 

Fig. 1. — uvby metallicities for planet-bearing stars obtained with our photometric cali- 
bration in comparison with spectroscopically determined metallicities from the CSR 2001 
compilation; solid line is a 1:1 correlation 

Fig. 2. — Distribution of residuals between photometrically and spectroscopically determined 
metallicities, with Gaussian fits to the distributions overplotted in the same line style. Our 
new photometric calibration is shown as a solid line; the SN calibration is shown as a dashed 
line. 

Fig. 3. — Distribution of residuals between photometrically and spectroscopically determined 
metallicities, with Gaussian fits to the distributions overplotted. Only measurements yielding 
a spectroscopic determination of [Fe/H] > 0.0 are plotted. Our new photometric calibration 
is shown as a solid line; the SN calibration is shown as a dashed line. 

Fig. 4. — Distribution of calculated values of planeticity for the 73 nearby Lick Survey stars. 
Stars without planets are shown as a solid line; stars with confirmed planets are shown as a 
filled histogram. 
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